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Lagrangian mechanics uses d'Alembert's principle of zero virtual work as an important starting 
point. The orthogonality of the force of constraint and virtual displacement is emphasized in liter- 
ature, without a clear warning that this is true usually for a single particle system. For a system 
of particles connected by constraints, it is shown, that the virtual work of the entire system is zero, 
even though the virtual displacements of the particles are not perpendicular to the respective con- 
straint forces. It is also demonstrated why d'Alembert's principle involves virtual work rather than 
the work done by constraint forces on allowed displacements. 
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I. INTRODUCTION 

The principle of zero work by constraint forces on virtual displacement, also known as d'Alembert's principle, is an 
important step in formulating and solving a mechanical problem with constraintsi*Siiii. In the simple systems widely 
used in literature, e.g., a single particle rolling down a frictionless incline, or a simple pendulum with inextensible 
ideal string, the force of constraint is perpendicular to the virtual displacement. This results in zero virtual work 
by constraint forces. It is often tempting to assume that the constraint forces are always orthogonal to respective 
virtual displacements, even for a system of particles. d'Alembert's principle then seems to be a consequence of this 
orthogonalitji^. 

In this article we study two simple systems: a double pendulum and an A^-pcndulum. In these systems it is observed 
that, the virtual displacements are not perpendicular to the respective constraint forces acting on individual particles 
(pendulum bobs). However, d'Alembert's principle of zero virtual work still holds for the systems as a whole. In these 
problems, the principle of zero virtual work is a consequence of, (i) the relation between the virtual displacements of 
coupled components (neighbouring bobs), and (ii) the appearance of (Newtonian) action- reaction pairs in forces of 
constraint between neighbouring particles. Thus d'Alembert's principle is more subtle and involved than it is often 
thought to be. Greenwood has rightly said, "... workless constraints do no work on the system as a whole in an 
arbitrary virtual displacement. Quite possibly, however, the workless constraint forces will do work on individual 
particles of the system"—. 

II. DOUBLE PENDULUM 

Let us first consider a double pendulum, with inextensible ideal strings of lengths L\ and L 2 - Let ri and v 2 denote 
the instantaneous positions of the pendulum bobs V\ and P 2j with respect to the point of suspension (see figure 1). 
When the system is suspended from a stationary support; the holonomic, scleronomous constraint equations are, 

M=Li, |n-r 2 |=i 2 . (1) 

The equations for allowed velocities are obtained by differentiating JQ) with respect to t, 

ri-vi=0, (r a -ri)-(v 2 -vi) = 0. (2) 

Thus, the allowed velocity vi of V\ is orthogonal to its position vector ri. The relative velocity of the second bob 
(V2 — vi) with respect to the first, is orthogonal to the relative position vector r 2 i = (r 2 — ri). Let rii and A2 be unit 
vectors orthogonal to (ri ) and (r 2 — ri ) respectively, 

ni • ri = and n 2 • (r 2 - ri) = 0. (3) 

From |J2J and we obtain the allowed velocities as, 

vi = &ini, v 2 — Vi = b 2 n 2 => v 2 = vi + b 2 n 2 (4) 
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where b\ and 6 2 are arbitrary real constants, denoting the magnitude of the relevant vectors. 




FIG. 1: Double pendulum: (a) position vectors, (b) orthogonality of tensions to part of virtual displacements, (c) cancellation 
of part of virtual work related to action-reaction pair (tension) 

The virtual velocities are defined as a difference between two allowed velocities, v k = v k — v^— . The allowed 
displacements dr k , and the virtual displacements Sr k are then given by, 

dri = vi dt = dqihi, 

5ri = dri — dr[ = (dqi — dq[)hi = Sqihi. (5) 
dr2 = V2dt = Vidt + & 2 n 2 <ii = dri + dg 2 n 2 , 

Sr 2 = {dri + dq 2 h 2 ) - (dr[ + +dq' 2 n2) = Sri + <5<7 2 n 2 . (6) 

It may be noted that, under the given holonomic, sclcronomous constraints (velocity and time independent), the sets 
of allowed and virtual displacements are equivalent. 

A set of allowed displacements {dri,dr 2 } is obtained by a specific choice of the numbers {61,62} or {dgi,dq 2 }. By 
making different choices of the set {dqi, deft} we get a whole family of allowed displacements, {dr^}, where, 

{dr k } {dri,dr 2 } {bl , b2} = {dr 1: dr 2 } {dqudq2 }. 

Or more precisely, 

{dr k } => {{dri,dr 2 } {blM} \bi,b 2 e 

where is the set of real numbers. Similarly, by choosing different set of quantities {5qi,Sq2} we get the family of 
virtual displacements, {<5rfc}, where, 

{Sr k } => {ST 1 ,5r 2 }{ SquS g a }. 

Thus, it is easy to see that for holonomic, scleronomous constraints, the set of all possible allowed displacements 
is the same as the set of all possible virtual displacements. This is in agreement with the fact that for holonomic, 
scleronomous systems, the sets {dr k } and {<5i"fc} satisfy the same equations^, namely, 

ri • dri = 0, (r 2 - r x ) • (dr 2 - dri) = 

ri • Sri = 0, (r 2 - ri) • {Sr 2 - Sri) = 

As the pendulums arc suspended by inextcnsiblc ideal strings, one may assume that the tensions in the strings act 
along their lengths. This essentially implies that there is no shear in the string to transmit transverse force. Thus the 
tension Ti is along (— ri) and tension T 2 is along (ri — r 2 ). Hence, the virtual displacement Sri for the first pendulum 
is perpendicular to the tension Ti, but the virtual displacement (5r 2 of the second pendulum is not perpendicular to 
T 2 . 



Ti • Sri = 

T 2 • Sr 2 = T 2 • (Sri + Sq 2 h 2 ) = T 2 • Sn + (T 2 • n 2 ),5 g2 = T 2 • 5n 



(7) 
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At this stage one may appreciate that Ti is not the entire force of constraint on V\. As a reaction to V\ pulling Vi 
with a tension T2, the second bob V 2 pulls the first bob V\ with a tension (— T2). Thus the virtual work by constraint 
forces acting on V\ and V 2 are given by, 

SWi = Ri • Sn = (T x - T a ) • Sr 1 = -T 2 • Sr 1 ^ 

SW 2 = R 2 ■ Sr 2 = T 2 • (5n + Sq 2 h 2 ) = T 2 • Sn ^ 0. (8) 

Although neither SWi nor SW 2 is zero, their sum adds up to zero. Thus the "equal and opposite" Newtonian reaction 
comes to our rescue, and we have a cancellation in the total virtual work. 

5Wi + SW 2 = Ri • 6ti + R 2 • 5r 2 = -T 2 • Sn + T 2 • Sn = (9) 

This shows that in the case of a double pendulum with stationary support, d'Alembert's principle utilizes the equal 
and opposite nature of the tensions between neighbouring bobs (Newtonian action-reaction pair). 

Let us now consider the double pendulum with a moving point of suspension. This gives us a system with a 
rhconomous constraint. Let the velocity of the point of suspension be vo- The constraint equations in this case are, 

|ri-v t|=Zi, \r 2 -n\=L 2 (10) 

The only non-trivial modification is that, the virtual displacements are no longer equivalent to the allowed displace- 
ments. 

For the first bob Vi, 

dri =vidt = vo dt + dqini, 

Sri = dvi - dr[ = (v\ c^ + dqin±) - (^g$ + dq[n{) = (dqi - dq[)&i = Sqini. (11) 

Therefore virtual displacement <5ri is a vector along ni, whereas allowed displacement dv\ is sum of a vector along 
ni and a vector along Vo . For the second bob V 2 , 

dr 2 = w 2 dt = Vidt + b 2 ii 2 dt = v Q dt + dqi&i + dq 2 h 2 

Sv 2 = (%&p + dqihi + dq 2 n 2 ) - (\\ c^ + dq[hi + dq' 2 h 2 ) = Sqihi + Sq 2 h 2 . (12) 

Thus dri and dr 2 arc not equivalent to Sri and 5r 2 . However the relation between ^ri and Sr 2 remains the same as 
in the case of a double pendulum with stationary support. 

Sr 2 = Sri + Sq 2 h 2 (13) 
Hence the above inferences, in particular, (0), ijHJ, © are true even in this case. 

III. iV-PENDULUM 

It is instructive to repeat the above exercise for a system of A^-pendulum joined end to end by inextensible, ideal 
strings. Let ri, r 2 , . . . , rjy denote the instantaneous position vectors of pendulum bobs Vi,V 2 t ■ ■ ■ , Tn respectively, as 
shown in figure [2 The constraint equations for this system are, 

|ri| = ii, |r fe -r fc _i| = L k , k = 2, 3, . . . , N. (14) 

The equations for allowed velocities, obtained by differentiating the above equations, are 

ri -vi=0, (r fc -r fc _i)-(v fc -v fc _i) = 0, k = 2, 3, . . . , N. (15) 

Let us introduce unit vectors ni, n 2 , . . . n^r, where is normal to the relative position of the bob Vk with respect to 
Vk-i, i-e., (r fe - r fe _i). 

n 1 -n=0, n fc -(r fc -r fc _ 1 ) = 0, k = 2, 3, . . . ,N. (16) 

From ()15|l and l|16(l . the allowed velocities are given by 

vi=&ifii, v fe - v fc -i = 6 fc n fe v fe = v fc _! + b k h k , k = 2, 3, . . . , N. (17) 
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(a) (b) (c) 



FIG. 2: JV-pendulum: (a) position vectors (position of the point of suspension in relation to the pendulum is schematic), (b) 
orthogonality of tensions to part of virtual displacements, (c) cancellation of part of virtual work related to action-reaction pair 
(tension) 

where {bk,k = 1, . . . ,N} are a set of real constants denoting the magnitude of the relevant vectors. The allowed 
displacements dr^ and the virtual displacements Srk are, 

dri = bihidt = dqi&i, Sri = Sqi&i, 

dr 2 = dri + b 2 n 2 dt = dr± + dq 2 n 2l 5r 2 = Sr ± + Sq 2 h 2 

drjy = drjy-i + fr/yfi/vcft = dr/v-i + dq^iiN, Sr^ = Sr^-i + Sq^iiN 



As is noted in the previous section for case of double pendulum, due to the holonomic, scleronomous nature of the 
constraints, the sets of allowed and virtual displacements are equivalent. From the above equations one can see that 
the virtual displacement of each pendulum (Srk) is a vector sum of the virtual displacement of the previous pendulum 
(5rk-i) and a component along the unit normal n^. Hence the virtual displacements (with the exception of <5ri) are 
not orthogonal to the corresponding relative position vectors. 

Let us now consider the constraint forces on each individual pendulum bob Vk- The bob Vk is pulled towards its 
point of suspension (the previous bob Vk-i) by a tension Tfc along (rfc_i — rfc). The next pendulum bob, Vk+i, is 
pulled towards Vk by a tension T^+i along (r^ — r^+i). In response to this, a reaction force (— T^+i) acts on the bob 
Vk along (rfc+i — r^). Thus between any two neighbouring pendulum bobs, there exists a pair of equal and opposite 
action-reaction forces. The total force on Vk is (Tfc — Tk+i) for k = 1, 2, . . . , N — 1. However, for the last pendulum 
Vn, the net constraint force is Tjy. 

The virtual work done by the constraint forces at different system points (particle positions) are, 

SWi = Ri • Sri = (Ti - T 2 ) • 5n = Ti • Sr x - T 2 ■ Sr ± 

5W k = R fc • Sr k = (T fc - T fc+1 ) • Sr k = T k ■ (Sr k -i + Sq k h k ) - T k +i ■ Sr k k = 2, . . . , N - 1 

SW N = Rjy • Sr N = T N ■ (Sr N -i + 5q N n N ) (18) 

As the strings of the pendulums are ideal, the tension Tfc acts along the length of the string, i.e., (rfc_i — rfc). Thus 



5 



the tension is normal to the unit vector n^. The above virtual work elements become, 

<fWi = -Ta-*ri, 6W k = T k -5r k - 1 -T k+1 -5r k , 5W N = T N ■ Sr N ^. (19) 

It is clear that the virtual work at each system point V k is non-zero. However if we sum the virtual work at all these 
system points, we observe a mutual cancellation and the total virtual work vanishes. 

Ef=i SW k = -T a • *n + (T 2 • 5ti - T 3 • Sr 2 ) + 

+ (Ti-i ■ 5r 3 -_ 2 - Tj ■ St^x) + (T, • St,-! - T j+l ■ Sr 3 ) + (T J+1 • 5r 3 - T j+2 ■ 5r j+1 ) + ... 

+ (Tjv-i • Svjf-a - T N • ^rjy-i) + Tjy • St N - X = (20) 

This vanishing of total virtual work is a consequence of (i) definition of virtual displacement, (ii) appearance of 
action-reaction pairs in the forces of constraint. The virtual work connected to each bob V k , is composed of three 
parts, (i) virtual work by the tension T k on the component of virtual displacement Sq k iik orthogonal to the relative 
position vector, (ii) virtual work by the tension on part of the virtual displacement 5r k -\ related to that of the 
previous bob, and (iii) virtual work by the reaction tension (— Tfc+i) (acting towards the next bob V k +i) on the 
virtual displacement 8r k — Sr k -i + 5q k iik- The first component for each V k vanishes because of orthogonality of the 
related force and virtual displacement. This is shown schematically in figure Efb). Due to the "equal and opposite" 
nature of action reaction pairs, and existence of a common term in the virtual displacement of neighbouring bobs, the 
other terms for each bob cancel with the related terms of its neighbours. Shaded areas in figure |3fc) illustrate this 
cancellation. 



TABLE I: Virtual work and d'Alembert's principle for simple and TV-pendulum 



System 

Pendulums with fixed string length 
holonomic constraints 


Stationary support 
scleronomous 
(Sr k ~ dr k y 


Moving support 
rheonomous 
(Sr k ^ dr k )* 


Simple pendulum 

R = T 


T ■ Sr = 


T ■ dr / 
R ■ dr 




T ■ Sr = 




R ■ Sr = 


R ■ Sr = 






JV-pendulum 

R& = (T fc - T k+1 ), k = 1, ...N - 1 
Rjv — Tjv 


T fc • Sr k jt 

N 

£ T fe -«5r fe /0 

fc=l 


T fc • dr k 

f; T fe ■ dv k jt o 

fe=i 


T fe ■ Sr k 

N 

£ T k -Sr k ^0 

k=l 


Rfc • Sr k 


Rfc ■ dv k 
J2 R fc -dr fe /0 

fe=I 




R fc ■ 5r k ± 




N 

J2 B, k -5r k =0 

k=l 




J2 R fe -5r fc =0 

fc=l 













' : Sr k and dr k are equivalent if the constraints are both holonomic and scleronomous. 

* : 5r k and dr k are not equivalent if the constraints are non-holonomic and/or rheonomous. 



Let us now study the iV-pendulum when its point of suspension is moving with velocity Vo- The system now has a 
rheonomous constraint as well, 

|ri-vot|=ii, \r k - r k -i\ = L k , k = 2, 3, . . . , N. (21) 

In accordance with the previous section, for /V-pendulum with moving support, the only non-trivial modification is 
that, the virtual displacements are no longer equivalent with the allowed displacements. For the first bob Vi, 

dr\ =vidt = vq dt + dq±TLi, 

Sri = dri - dr[ = {^p&p + dqin±) - {\>(tp + dq[iii) = {dqi - dq[)ni = Sqini. (22) 
Thus 8ti is a vector along rki, whereas dri is sum of a vector along ni and a vector along Vq. For subsequent bobs 



6 



dr k = v k dt = v k _ 1 dt + b k ii k dt = v dt + (%fii + dq 2 h 2 + h b k h k dt 

Sr k = {^)df. + dqxiix + dq 2 h 2 H h dq k h k ) - (^df + dq[hi + dq' 2 n 2 H h dq' k h k ) 

= <5<?ini + 5q 2 h 2 H h Sq k h k . (23) 

Thus 6r k and dr k are not necessarily equivalent. However the relation between 5r k -i and Sr k remains the same as in 
the previous case of N pendulum with stationary support. 

Sr k = (5r fc _i + Sq k h k (24) 

Hence the above inferences, in particular, (|18|) . 1 (19(1 . 120(1 are true even for an A-pcndulum with moving support. 
The adjacent table summarizes the results presented in sections II and III. 



IV. CONCLUSION 



The zero virtual work principle of d'Alembert identifies a special class of constraints, which is available in nature, and 
is solvable^. Two noteworthy features of d'Alcmbcrt's principle are, (i) it involves the virtual work (R k -Sr k ), i.e., work 
done by constraint forces on virtual displacement 5r k and not on allowed displacement dr k , and (ii) the total virtual 
work for the entire system vanishes, i.e., ^(Rfc-iJrfc) = 0, though virtual work on individual particles of the system need 
not be zero (R k ■ 5r k ) ^ 0. For holonomic (velocity independent) and sclcronomous (time independent) constraints, 
e.g., pendulum with stationary support, the allowed and virtual displacements are collinear and hence a distinction 
between work done on allowed displacement and that on virtual displacement is not possible. For understanding the 
nature of this distinction one needs to study a system which is either non-holonomic or rheonomous. A pendulum with 
moving support, or a particle sliding down a moving frictionlcss inclined plane arc examples of simple rheonomous 
systems. In constrained systems involving a single particle, the second feature mentioned above, in reference to 
d'Alembert's principle, becomes irrelevant. As there is only one particle, there is no summation in virtual work, and 
(R • Sr) = 0. This implies that the force of constraint R is normal to the virtual displacement Sr. In order to really 
appreciate the importance of summation in the total virtual work, one needs to study system of particles involving 
several constraints. The double pendulum and iV-pendulum, particularly with moving support, present two of the 
simplest systems illustrating the subtlety of d'Alembert's principle. 
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